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We show that in the severe slowing down temperature regime the relaxation of antiferromag- 
netic rings and similar magnetic nanoclusters is governed by the quasi-continuum portion of their 
quadrupolar fluctuation spectrum and not by the lowest excitation lines. This is at the heart of the 
intriguing near-universal power-law temperature dependence of the electronic correlation frequency 
Uc with an exponent close to 4. The onset of this behavior is defined by an energy scale which is 
fixed by the lowest spin gap Aq. This explains why experimental curves of uJc for different cluster 
sizes and spins nearly coincide when T is rescaled by Ao. 
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I. INTRODUCTION 

A central issue for the control and manipulation of elec- 
tronic spin degrees of freedom in the field of molecular 
nanomagnetism^ and related areas is the understanding 
and characterization of the various microscopic mecha- 
nisms of relaxation and decohcrencc which stem from the 
interactions with the underlying degrees of freedom of the 
host lattice. It is by now well established that in the large 
majority of magnetic clusters, the relaxation of the mag- 
netization shows a dramatic slowing down which sets in 
already at relatively high temperatures T and is char- 
acterized by a single electronic frequency cutoff ujc^'"^'^. 
While the majority of freezing mechanisms reported in 
the literature [e.g., onsite anisotropy in single molecule 
magnets^, anisotropy and critical fluctuations in single 
chain magnets^, phonon trapping in Niio^], have been 
understood to a large extent on the basis of the nature 
of the lowest excitations, the case of antiferromagnetic 
rings (AFMR's) has been exceptional and very intrigu- 
ing. Indeed, as shown by Nuclear Magnetic Resonance 
(NMR) experiments by Back et alJ' , in the regime where 
the dramatic slowing down effect takes place, lJc follows a 
near-universal power-law T dependence with an exponent 
close to 4, a behavior clearly inconsistent with a relax- 
ation scenario based on the lowest excitation lines. It is 
also striking that the experimental curves of lOc from var- 
ious magnetic rings with different sizes and spins s > 1/2 
nearly coincide when T is rescaled by the lowest spin gap 
Ao. 

On the theoretical side, the first microscopic calcula- 
tion of ijJc in AFMR's has been given by Santini et al."^ 
who conclude on this problem that the relevant T regime 
is too narrow to substantiate any power-law scaling and 
that instead lOc is approximately exponential in 1/T. 

Here we present a new microscopic theory which sup- 
ports the scaling hypothesis of Back et al.^ and resolves 
its physical origin in a clear and transparent way. Our 
central key finding is that, in the severe slowing down T 



regime, the relaxation is not governed by the lowest exci- 
tation lines of the quadrupolar spectral density but by a 
forest of quasi-continuum excitations at higher energies 
(cf. Fig. 2). This is at the heart of the power-law temper- 
ature dependence of coc with an exponent which is close 
to 4 (cf. Eq. (12)). The onset of this behavior is defined 
by an energy scale which is fixed by the lowest spin gap 
Aq , and this explains why curves from different rings fall 
almost on top of each other when T is rescaled by Ag. 
The present theory is corroborated by a model calcula- 
tion of the nuclear spin-lattice relaxation rate 1/Ti which 
includes one-phonon, two-phonon as well as Raman pro- 
cesses and shows excellent agreement with experimental 
data for the Crs, FceLi and FegNa clusters. 

Our theory is based on the formalism developed in Ref. 
[3] . The advantage of this method is that it is based on an 
analytical formula which gives lOc in terms of a frequency 
overlap between an electronic and a phononic spectral 
density function [cf. Eq. (3) or (5) below]. Thus by mon- 
itoring the variation of these spectra with T and lo one is 
able to identify the dominant relaxation channels in each 
T regime of interest. The central approximation for the 
derivation of this analytical formula of is that the spin 
Hamiltonian TIq commutes with the total magnetization 
Sz of the cluster. For the majority of magnetic clusters 
reported in the literature this approximation is an excel- 
lent one for the study of relaxation phenomena at not too 
low T (T > 1 K). Indeed, the dominant energy scale in 
the problem is set by the isotropic Heisenberg exchange 
and very often by a uniaxial onsite anisotropy. There 
exist several types of anisotropic interactions which do 
not conserve Sz, but these are typically very small (< 1 
K) and thus play a significant role only at very low T. 
Thus, although the present work is devoted to the slow- 
ing down mechanism in AFMR's the method presented 
here can be applied to the majority of magnetic clusters 
reported in the literature. In particular, we will argue 
that the slowing down mechanism in AFMR's must be 
common in the general class of antiferromagnetic nano- 



2 



magnets (which includes e.g., grids and other clusters) 
due to their very similar spectral structure (cf. below). 



II. METHOD AND MODEL 

We consider a magnetic ring with an even number N 
of spins s > 1 described by the Hamiltonian 



Ho = J^Si ■ Sj+i + gusBS, 



(1) 



with periodic boundary conditions. The first term de- 
scribes the antiferromagnetic (J > 0) exchange between 
nearest-neighbor spins, the second is the Zeeman energy 
in a field B ~ Be^, S = J^i is the total spin, g ~ 2, 
and /iB is the Bohr magneton. It is known'^'®'^ that the 
energy spectrum of Tio is bounded from below by an ex- 
citation band Es ~ Ao5(5 + l)/2, where Aq ~ AJ/N 
(cf. Rcf. [9] for the validity range of this scaling) is the 
lowest singlct-triplct gap. At higher energies there ap- 
pear a forest of quasi-continuum excitations which set in 
progressively above the lowest band. This dense spectral 
structure above a certain energy scale is very common in 
finite unfrustrated antiferromagnets^°'^^ and, as we show 
below, affects the relaxational behavior in a very charac- 
teristic way. 

We are interested in the damping of the equilibrium 
fluctuations of the total magnetization Sz at not too low 
T and for huje = gfJ-sB <^ J. This damping is trig- 
gered by the (phonon-drivcn) fluctuating portion of vari- 
ous anisotropics. Here we consider the quadrupolar spin- 
phonon channel^'^^'^^'^'' which we write in the general 
form 



N 



N 



's-ph = 5]Q(s0-*(r.0 = E E 



(2) 



emission coefficient J, 



(w') of the host lattice 
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As we show below, this formulation is very fruitful since 
it allows to identify the frequency regime (or the excita- 
tions) which gives the dominant contribution to lUc in the 
severe slowing down T regime. 

To proceed we shall make the reasonable assumption 
that the strain fields (e.g., local librations of the lig- 
and groups) are uncorrelated between different magnetic 
sites, i.e., 



J. 



(4) 



(below we drop the site indices since all operators shall 
refer to a single site). We further note that the SU(2) 
invariance of TIq at B = necessitates that Jg gt (uj) = 

jQ_^Q-t^{io)dmm' ■ Since this remains true at finite B for 
ksT > huje, we may replace Eq. (3) by 



dco' 4^J^t^(c.') , (5) 



where Js(^) — Jn n^'i^) 



m 



(6) 



which defines implicitly an average coupling field <i>(r,t). 

For our purposes we use a simple Debye model consist- 
ing of 3 acoustic branches with a common sound velocity 
c and a Debye cutoff lud (cf. below), and keep from the 
strain tensor^^'^'^'^^ e only its isotropic (scalar) portion 



r/(r, i) = V 



2Mc 



kcr 



gi(kT+Wk„t)(jt 



h.c. 



(7) 



where $(ri) arc functions of the local strains or rota- 
tion fields, and Qmi^i) are the quadrupolar operators 
Q±2(s) = s±, Q±i(s) = {s±Sz+SzS±)/2, and Qo{s) = s\. 

It has been found, both experimentally^'^ and 
numerically^, that the damping of Sz is mono- 
exponential (or Markovian) in a large number of nano- 
magnets. The physical origin of this central feature has 
been shown'^ to arise from a dynamical decoupling of Sz 
from the remaining slow degrees of freedom which, in 
turn, follows from the discreteness of the energy spec- 
trum and the conservation law [52,7^0] = 0. According 
to the general expression (27) of Ref. [3], lOc is given by 



where ihFi, 



\Sz 1 Qi 



ml 
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, /3 



(3) 

l/fc^T, and 



Xo = Pi^Sz) is the isothermal susceptibility. Thus ujc is 
proportional to the frequency overlap between the ab- 
sorption {-(jj' < 0) coefficient J^. (— "^O and the 
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(a) One-phonon (b) Two-phonon (c) Raman 

FIG. 1: (Color Online) Schematic representation of the three 
lowest order processes which contribute to coc for any given ab- 
sorption frequency lj of the nanomagnet. (a) The one-phonon 
absorption process is proportional to the average phonon en- 
ergy density cjp(cj)n(tj), where p(w) is the phonon density of 
states and n(uj) = (exp{Phw) — 1)~^ is the Bose-Eistein dis- 
tribution function, (b) The two-phonon absorption process is 
proportional to the integral / duj' uj' p{uj')n{uj'){uj — uj')p{ijj — 
Lj')n(uj — ui'). (c) The inelastic scattering Raman process is 
proportional to the integral J duj' lo' p{uj')n{ijj')(u)' — uj)p{uj' — 
Lo)[n{w' -(jo) + 1]. 
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where u(r, t) denotes the displacement field, M is the to- 
tal mass of the crystal, cr runs over the three polarization 
states and tJkcr = c|k|. Wc may then expand $(r, t) as 

$(r,t) ~ wi7?(r,t) +W2?7^(r,t) , (8) 

where vi and V2 define two spin-phonon coupling energy 
parameters whose importance will become clear below. 
Since Js(w) is sharply peaked at the Bohr frequencies 
UJB of the cluster, i.e., Jsicu) = J^^^ Jsi'^B)S{uj - ujb), 
we may rewrite Eq. (5) as 

(9) 

Using the statistical factor n{u)) = {e^^ — 1) ^ and the 
mass density pm, we have for < a; < wd: 

3 ^ 

Jrmi^) = ^n{uj)uj^ = Jip/i(tj) (10) 

which is the contribution from direct onc-phonon pro- 
cesses, while J^2jj2(lo) is the sum of the two-phonon and 
Raman contribution, given respectively by 

J2ph{to) = ^ dijj'oj'^n{uj')[LO - uj')'^n{ijj - uo') 
Jo 

Jr{uj) = 2A duj'Lu'^n{uj') x 

J UJ 

(lu' ~ Lof [niu' - LO) + 1] , (11) 

where A = ^fg-p^prr- These processes are represented 
schematically in Fig. 1. 

Equations (9-11) arc the starting point of our calcula- 
tions for LOc- The quadrupolar spectral density Js(w) is 
obtained by a full thermodynamic calculation using exact 
diagonalizations . 



III. SLOWING DOWN EFFECT AND 
NEAR-UNIVERSAL POWER-LAW SCALING 

Figure 2 shows Js(-tj), Jiph(w), J2ph{uj) and Jr{uj) 
for Crg at hjJe = O.O66A0 and for fc^T = 0.19Ao (a) and 
2.22A0 (b). The values of the parameters used here are 
the ones that fit the l/Ti data (cf. Table I). Apart from 
identifying a number of excitation lines such as Li (tran- 
sitions within the lowest triplet), L2 (transitions from the 
lowest triplet to the lowest quintet) and L3 (transitions 
from the lowest singlet to the lowest quintet) we find 
a number of general features which originate in the over- 
all spectral structure of AFMR's described above. We 
first emphasize the increasingly gapped structure at low 
Lo and T, and the fact that Js(— w) has essentially no 
weight below the line L2. This structure is altered very 
quickly by thermal excitations since the latter increase 
the number of available resonant channels and thus give 



rise to a dense fluctuation spectrum at higher T . Impor- 
tantly, Js(— w) remains appreciable over an overall band- 
width huJraa.^ ~ 10 — 20Ao. Givcu uow the behavior of 
the phononic density (cf. Fig. 2) the relaxation process 
can be understood as follows. At fc^T < Tiw,, ^ Aq the 
system is essentially opaque to the available (thermally 
excited) phonons since there is no appreciable overlap 
between Js(— w) and Jiphiuj). Only the Zecman line Li 
contributes, giving Wc ^ ujl/{e^'^'' — 1). However, as 
soon as the one-phonon spectral peak reaches the quasi- 
continuum regime (above the line L2) a significant num- 
ber of spin-phonon resonant channels are quickly ther- 
mally activated. In fact. Fig. 2(b) shows that Jiphi}-^) 
samples almost the entire electronic spectral weight up 
to Wmax already at fc^T ~ 2Ao. This marks the exis- 
tence of a special regime 0.2Ao < fcsT < 2Ao, where a 
dramatic change of by 3-4 decades (cf. Fig. 3) takes 
place. We emphasize here that the relaxation process in 
this T regime is driven by the quasi-continuum portion 




FIG. 2: (Color online) The origin of the dramatic slowing 
down effect in AFMR's, demonstrated here for the Crg cluster 
at ksT = O.I9A0 (a) and 2.22Ao (b), with fiuj^ = O.O66A0 
and hujD ~ lOJ. The correlation frequency Wc is proportional 
[cf. Eq. (3) or (5)] to the overlap between the quadrupolar 
spin density Js{—u}) (series of (5-peaks, here shown with a 
logarithmic mesh in tj) and the phononic density J-^f-^{i^) 
which includes one-phonon (solid red), two-phonon (dotted 
red) and Raman (dashed red) processes. We also show (solid 
blue line) the integrated density f^duj'Js(—Lj'). The lowest 
three groups of excitation lines at low T are denoted by Li, 
L2 and L3 (cf. text). 




FIG. 3: (Color Online) Near-universal behavior of ujc vs. 
fcflT/Ao for rings with different N and s. The dramatic 
slowing-down (by 3-4 decades) takes place in the regime 
0.2Ao < fcsT < 2Ao, where uJc follows Eq. (12). The lim- 
its of this distinctive T regime can be understood by looking 
at Fig. 2: It begins when the one-phonon spectral peak starts 
overlapping the quasi-continuum portion of the quadrupolar 
density, and terminates when it has sampled its entire band- 
width. At fcsT > 2Ao, LJc oc T (Raman processes or high- 
frequency phonons alter this behavior), while at very low T, 
LOc ~ cJeAe''*"^ - 1) (cf. text). 
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FIG. 4: (Color Online) Comparison between theory (solid red 
lines) and experimental data (symbols) for R [cf. Eq. (14)] 
and u)c for Crg (left) and FeeLi (right). The one-phonon 
(dashed), two-phonon (dotted) and Raman (dash-dot) con- 
tributions to LOc are also shown, along with the respective 
power-law scalings (dash-double dot, blue lines) for 0.2Ao < 
ksT < 2Ao. 



of the quadrupolar spectrum and is minimally affected 
by the lowest excitation lines. This shows that lUc is 
not a nearly exponential function in 1/T in this regime. 
In fact, we may go one step further and inquire into 
the observed power-law scaling by employing a steepest- 
descent expansion in Eq. (5) which relies on two central 
ingredients: (i) the dominant contribution to the over- 
lap comes from the regime around the one-phonon peak 
Hujp{T) ~ 2.82kBT (note that Fig. 2 is in logarithmic 
scale), and (ii) the quasi-continuum character of Js{—lu) 
in the respective T and lo regime. To leading order 

LOc (X rV(T) (12) 

where /(T) = Js{-oOp{T),T)/xoT, and the strong 
pre-factor emerges from the functional form of Jiph{^, T) 
alone. The function f{T) varies with N and s but shows 
always a much weaker (sub-linear) T dependence. This 
can be seen by our calculations shown in Fig. 4 (cf. text 
below) which are in excellent agreement with experimen- 
tal data. Hence, the y3-5±o.5 gcaling law reported in Ref. 
[5] is largely due to the above leading pre-factor. It is 
thus essentially a characteristic fingerprint of one-phonon 
acoustic processes and of the quasi-continuum nature of 
the electronic quadrupolar spectrum at high energies. 

All of the above features are demonstrated in Fig. 3 
which shows the one-phonon contribution to lOc for vari- 
ous N and s. Indeed, one identifies a distinctive regime 
0.2Ao ^ ksT < 2Ao, where lOc shows a dramatic drop by 
3-4 decades. Figure 3 explains yet another of the central 



findings of Baek et al.^, namely that all curves fall almost 
on top of each other when plotted against fesT/Ag. This 
is clearly due to the overall similar spectral structure of 
AFMR's and, in particular, due to the fact that the for- 
est of quasi-continuum excitations sets in at an energy 
scale fixed by Aq ~ AJ/N (and not by J). 

We should add here that the details of the low-energy 
spectrum such as the character of the lowest L band or 
the spin-wave E band^ do not play any special role com- 
pared to other excitations. The quasi-continuum portion 
of the spectrum contains excitations from and towards 
both the L and E bands, but it is the global dense aspect 
of the spectrum that matters. 



IV. COMPARISON TO NUCLEAR 
SPIN-LATTICE RELAXATION RATE DATA 

Let us now describe our calculations for 1 /Ti and com- 
pare to experimental data for FeeLi, FegNa and Crg. We 
begin with the expression^'^'^* 

l/Ti = A,,Js„,soA^l) = '^Js^sA^l) (13) 

where Sq is the ionic spin with the shortest distance tq 
from the nuclear spin, and Azz = ^llln/i^t ^^'^ '^o^" 
responding hyperfine amplitude^^. Using'^ Js^sA'^l) = 
'^(^^z)z;^' and Xmoi = NAg^filPiSSl) (where Na is 
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TABLE I: Known data^"'^^ (first three columns) and fitting 
parameters for FegLi, FeeNa, and Crg. Here c = "1 5°^ . 



Cluster 


Pm(g/cm^) 




J(K) 




«2/c5(K) 


ro(A) 


Fee Li 


1.45 


217.8 


21 


0.498 


2.861 


4.59 


FeeNa 


1.42 


209.8 


28 


0.200 


1.277 


4.16 


Crg 


1.08 


154±10 


17.2 


1.124 


14.376 


4.10 



Avogadro's number) we get 



R: 



1 



TlXmolT 



4585.3 



(14) 



where l/Ti is given in ms, 7„ in MHz/T, Xmoi^ in cm'' 
K/mol, ro in A, and lUc.l in MHz. According to Eq. (14), 
the fit of ro is controUed by the magnitude of 1/Ti, while 
the values of Vi and V2 of Eq. (9) can be found by ad- 
justing the position and the width of the l/Ti peak. We 
should note however that V2 affects lJc only at the high-T 
side of the peak where high energy phonons start to play 
a role and thus its estimate is generally less accurate. 

Figure 4 shows our optimal fits to the data of R and lUc 
for Crg and FcgLi, along with the separate contributions 
of direct, two-phonon and Raman processes. Similar fits 
(not shown here) are obtained for FegNa. The corre- 
sponding estimates for ro, vi and V2 are given in Table 1 
and are of the correct order of magnitude. As to the value 
of the Debye cutoff, we find that our fits remain very good 
for fiior) ^ lOJ^^, which is consistent with the reported 
estimates of the Debye temperature 8d [cf. Table 1]. 
The agreement at intermediate T where the enhancement 
of l/Ti takes place is remarkably good. In particular, 
we find that the relaxation is dominated by one-phonon 
processes up to ksT ~ 2Ao for Crg {Ao/ks ^ 9.607K) 
and up to ~ 6A0 for FegLi (Aq/Zcb ~ 14.526K). On the 
other hand, the power-law scaling (dashed-double dot 
blue lines) is valid up to ~ 2Ao for both clusters. We 
also find that our model docs not account for the behav- 
ior at very high T. This can be either due to another 
contribution to l/Ti (e.g., electronic T2 processes"'^*'^'''^^ 
not included in Eq. (14)) or due to the actual details of 



the high-energy phonon modes which should play a role 
at high T. Finally, the discrepancy at very low-T in Crg 
arises from an additional peak in 1 /Ti which is currently 
not understood. 



V. CONCLUSIONS 

We have presented a microscopic theory which identi- 
fies the mechanism responsible for the dramatic slowing 
down effect observed in antiferromagnetic wheels. Our 
central key result is that in this class of nanomagnets the 
relaxation is driven by the quasi-continuum portion of 
the quadrupolar spectrum and not by the low-lying exci- 
tations. This is at the heart of the intriguing power-law 
T dependence of the electronic correlation frequency lOc 
with an exponent close to 4. The onset of this scaling 
is fixed by the lowest spin gap Aq ~ AJ /N (and not by 
J) and this explains why experimental curves of lJc for 
different cluster sizes and spins nearly coincide when T is 
rescaled by Aq. Since the spectral structure of AFMR's is 
very common in finite-size antiferromagnets most of the 
above qualitative features must carry over to other anti- 
ferromagnetic clusters as well. Hence we believe that the 
present slowing down mechanism of AFMR's is common 
in the more general class of antiferromagnetic clusters. 
More generally, the present theory can be applied to the 
majority of clusters reported in the literature for T > 1 
K. It can also be extended to clusters with s = 1/2 by 
using an appropriate relaxation channel.^ It is thus our 
hope that it will be valuable for further investigations 
on the understanding and characterization of relaxation 
mechanisms in magnetic nanoclusters. 
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